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SUMMARY 

A uniaxial approach was developed for calculating cyclic creep and stress 
relaxation at the critical location of a structure subjected to cyclic thermo- 
mechanical loading. This approach was Incorporated Into a simplified analyti- 
cal procedure for predicting the stress-strain history at a crack Initiation 
site for life prediction purposes. An elastic finite-element solution for the 
problem was used as Input for the simplified procedure. The creep analysis 
Includes a self-adaptive time Incrementing scheme. Cumulative creep Is the 
sum of the Initial creep, the recovery from stress relaxation and the Incre- 
mental creep. The simplified analysis was exercised for four cases Involving 
a benchmark notched plate problem. Comparisons were made with elastic-plastic- 
creep solutions for these cases using the MARC nonlinear finite-element com- 
puter code. 


INTRODUCTION 

This paper discusses the development of a uniaxial creep procedure for 
computing cyclic creep and stress relaxation at the critical location In a 
structure. While nonlinear, three-dimensional, finite-element computer codes 
are available for calculating cyclic stress-strain response, these are fre- 
quently too costly In computing time and other resources to be practical In 
the early design stages for complex structures. Computing costs In particular, 
limit three-dimensional, nonlinear finite-element analyses to a small number 
of cycles which may not be sufficient for life prediction purposes In cases 
Involving creep ratchettlng. The NASA Lewis Research Center has been sponsor- 
ing contractual and Inhouse efforts to develop simplified methods for predict- 
ing the stress-strain history at a crack Initiation site (ref. 1) as part of 
Its life Improvement programs for hot gas path components In aircraft engines 
and space power vehicles. This study was undertaken to Improve the creep 
analysis capability of an In-house developed simplified Inelastic analysis 
procedure. This procedure (refs. 2 to 3) performs a uniaxial nonlinear analy- 
sis for the critical location In a structure based on elastic solutions from 
three-dimensional finite-element analyses. 

The underlying assumption of this simplified procedure Is that the Ine- 
lastic regions are localized and constrained by the bulk of the surrounding 
elastic material. Therefore, the total strain history calculated from three- 
dimensional, elastic finite-element analyses for various points In the mission 


I 



cycle can be used as Input In elastic-plastic analyses for cyclic thermal 
loading. The elastic strain cycle Is divided Into a sufficient number of load- 
time steps to realistically determine the cyclic stress-strain hysteresis loops 
from Inelastic analyses. Neuber-type corrections (ref. 4), Including correc- 
tions for residual stresses, can be Incorporated Into the plasticity computa- 
tions to account for strain redistribution during mechanical loading cycles 
(ref. 5). This analytical procedure, with appropriate constitutive models to 
define the material cyclic behavior, was automated In a computer code. Since 
the simplified procedure Is uniaxial, the Input and output Is based on von 
Mlses effective stresses and strains. Signs have to be assigned to the Input 
effective strains from the elastic analyses by the user, usually on the basis 
of the sign of the dominant principal strain. A kinematic hardening model was 
used to determine cyclic yielding for the problems discussed herein. 

A power law type of creep constitutive equation was assumed to represent 
primary stage creep. This was used In conjunction with a time hardening rule 
for relating creep strain change to stress change. There Is also provision In 
the computer code for a modified strain hardening rule where the primary creep 
Is reinitiated and reversed with each reversal of the creep straining direc- 
tion. The unlaxlallty of the simplified procedure requires that the user 
exercise judgement In directing the form of the creep analysis. Three options 
are available In the code: (1) stress relaxation at constant strain, (2.) cumu- 

lative creep at constant stress or (3) a combination of stress relaxation and 
creep accumulation. This last option Is the one of Interest for most practical 
problems Involving high-temperature creep. Previously, the combined stress 
relaxation and creep strain change was obtained by averaging the results for 
the first two options for each specified time step. 

The purpose of this study was to develop a more rational and accurate 
creep analysis approach for combined creep and stress relaxation for Inclusion 
In the simplified analysis procedure. The approach that was developed Included 
a self-adaptive time Incrementing scheme for subdividing the specified dwell 
times Into Increments. At the start of the creep computations for each dwell 
time, the time period of the Initial Increment was automatically determined 
from preset tolerances on the permitted stress and Incremental creep strain 
changes. Creep for any subsequent Increment was calculated using a time period 
that was a constant multiple of the previous time period. In this approach, 
the cummulative creep strain at the end of the time period Is the sum of the 
creep at the beginning of the time period, the creep recovery from stress 
relaxation and the Incremental creep strain during the time period. The creep 
recovery component Is obtained by first calculating the stress relaxation for 
each Increment on the basis of constant strain boundary conditions, then assum- 
ing some proportionate recovery of the stress relaxation and calculating the 
creep strain change required for this recovery. For the Initial Increment, 
the creep recovery portion of the stress Is assumed to be two thirds of the 
stress relaxation. For each subsequent Increment, the percent of creep recov- 
ery Increases by a fixed amount until full recovery Is attained. These assump- 
tions were derived from experience In conducting a large number of nonlinear 
finite-element analyses of elastic-plastic-creep problems. 

The Improved simplified creep analysis was Incorporated Into the previ- 
ously developed computer code. Input consists of Inelastic material proper- 
ties In the form of coefficients and exponents for elastic-plastic and creep 
constitutive equations, dwell times at precrlbed points during the cycle and 
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the total strain history at the critical location from elastic finite-element 
analyses. Four cases Involving a benchmark notched plate problem (ref. 6) 
subjected to cyclic mechanical loading were analytically examined to verify 
the accuracy of the simplified creep analysis approach. The first two cases 
had maximum creep strains of about 0.2 percent during the tensile or loading 
part of the cycle and the third case had a maximum creep strain under 
0.1 percent during the compressive unloading part of the cycle. The fourth 
case Involved a combination of both tensile and compressive creep as an Illu- 
stration of a cyclic creep problem. Elastic and nonlinear finite element 
analyses were performed for all four cases using the MARC computer code. The 
elastic solutions for the notch root location were used as Input for the sim- 
plified analysis computer code. Comparisons between the simplified and MARC 
nonlinear finite-element analyses of these cases were made with particular 
attention to the creep portions of the stress-strain histories at the critical 
location. 

The comparisons demonstrated that the uniaxial creep analysis method can 
duplicate the creep strains In the cyclic stress-strain hysteresis loops com- 
puted from the nonlinear, three-dimensional, finite-element analyses to a high 
degree of accuracy and have greater computational stability than the latter. 
For the benchmark problems, the simplified procedure used less than 1 percent 
of the central processor (CPU) time required by the MARC cyclic analyses. 


SYMBOL LIST 

A,B,C temperature-dependent constants In creep power law, equation (1) 
E modulus of elasticity 

E e modified modulus of elasticity, 3 E/( 2(1 + p)) 

E p work hardening slope (fig. 3) 

K,n temperature- dependent constants In cyclic stress-strain law, 
equation (13) 

r time multiple constant 

t creep time 

t^ creep time at beginning of 1th time Increment 

a effective stress change tolerance coefficient 

time scale coefficient for 1th time Increment 
tc c \ equivalent creep strain change at 1th time Increment 
At^ 1th time Increment 

e equivalent total strain 
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e c equivalent creep strain 

c c -j equivalent creep strain at beginning of 1th time Increment 

e C R^ equivalent creep recovery strain at beginning of 1th time Increment 

e p equivalent plastic strain 

cp maximum equivalent plastic strain In cycle (fig. 3) 
percent creep recovery for 1th time Increment 
u Poisson's coefficient 

a effective stress 

a' maximum effective stress In cycle (fig. 3) 

effective stress at beginning of 1th time Increment 
oyi Initial yield stress In loading part of cycle 

o v 1 ' Initial yield stress In unloading part of cycle 


ANALYTICAL PROCEDURE 

A simplified cyclic creep procedure was developed for computing cyclic 
creep and stress relaxation at the critical crack Initiation critical location 
in a structure. In this approach, a power law type of creep constitutive 
equation was assumed to represent uniaxial primary creep behavior. 



( 1 ) 


A typical stress-strain cycle with creep dwell times at the maximum 
and minimum total strain parts of the cycle Is Illustrated In figure 1. A 
modified time-hardening law, as shown In figure 2, was proposed to determine 
the stress change at the 1th time Increment during a dwell time In the cycle. 
The equation, which Is derived In Appendix A, takes the form 
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where yi t the time scale coefficient, was assumed to Increase proportion- 
ately between the two limiting values (yi and y®) • These values were 
computed from the following expression 
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In equation (3), a Is the tolerance on the allowable stress change per 
Increment and r Is a constant relating the current time Increment to the 
preceding one. The latter Is defined as 


At 
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The value of r Is assumed to be 1.5, the same ratio used In MARC auto- 
matic creep computations. It follows that 


At^ +i = 1.5 At^ 


(5) 


and 


* 1+1 = *1 + At 1 


Equation (3) yields two limiting values: 
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To better fit data points with those obtained from MARC creep analyses, 
Instead of using equation (3) for yi . Y} was determined by a constant 
multiple of 7/6 of yi_i until Its value reached y«>- Then, for the 
remaining time Increments y» will be used throughout. 

In equation (2), denotes the percent of creep recovery at the 

1 th time Increment. For the Initial time period, the creep strain recovered 
Is assumed to be two-thirds of the stress relaxation. For each subsequent 
Increment, the percentage Increases by a fixed amount until a full recovery Is 
attained. Thus 

2 1 - 1 
6 1 = 3 * 30 


1=1 . . . . n 


(7) 


Full recovery Is assumed to have been attained at the tenth Increment. 

The above assumptions were derived from experience In conducting a large number 
of nonlinear finite-element analyses of elastic-plastic-creep problems. 


It should be noted that. In using equations (2) and (5) for each dwell 
time In a cycle, the Initial time period. At-], must be estimated first. A 
self-adaptive time-incrementing scheme (eq. (5)) was used to subdivide the 
remainder of the specified dwell time Into a number of Increments. At the 
start of the creep computation, the Initial stress level, a- 1 , was known from 
he calculated stress for the preceding Increment. The Initial time period, 

At-] , was estimated from the preset tolerances on the permlssable stress and 
Incremental creep strain changes, the Initial stress level, the Initial creep 
recovery, and the constants of the creep constitutive power law. The resulting 
equation takes the form 


At, 
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where Ae c Is the change of the Initial creep, which Is assumed to be equal 
to the change of the corresponding elastic strains. Since the von Mlses yield 
criterion was used, the total strain from a uniaxial stress-strain curve has 
to be converted Into a modified equivalent total strain (ref. 8). This modi- 
fied elastic part of the equivalent total strain corresponds to the measured 
elastic strain multiplied by 2(1 + yj/3- For computational convenience, this 
factor was Included In a modified elastic modulus defined by E e = 3E/(2(1 + p)) 
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Using the relationship (o-) - 02 )/E e for Ae c antl Imposing the stress change 
tolerance, equation (8) Is reduced to 



With at-] , known from equation (9), the subsequent time Increments follow 
equation (5). At the beginning of the dwell time; 1 =1, t] =0, e C ] =0 and 
at] was obtained from equation (9). At any other creep time, the cumulative 
creep strain up to t-j+i Is the sum of the creep strain, e c ^, at the begin- 
ning of time period t^, the creep recovery, QiecRI. f rom stress relax- 
ation, and the Incremental creep strain, a c c ^. These relations can be expressed 
as 

C c(1+1) = e c1 + Ac c1 " VcRI (10) 

where 

“cl {it A (t Hl - *?) <"> 

and 

‘cRt ' (if <Vl - ‘X <’ 2 > 

It should be noted that equation (12) Is exact only If the fully relaxed 
stress under constant strain boundary conditions Is used for rather 

than the partially recovered stress shown In figure 2. However, the difference 
between these two stress values Is slight for the small tolerance normally 
held on the stress change. When full recovery Is attained and 0^=1, these 
stress values become Identical. In view of these considerations and for the 
sake of simplicity, the partially recovered stress level shown In 

figure 2 may be accepted as the one to use In equation (2). 

The computational scheme started at the known effective stress-equivalent 
strain at the beginning of creep. The Initial time period was determined from 
equation (9) and the Initial time scale yi from equation (6). Subsequent 
time Increments follow equation (5). With these quantities known from the 
above procedure, the stress level at the end of the Initial time period can be 
obtained from equation (2) and the corresponding cumulative creep strains from 
equations (10) to (12). Subsequent stress relaxations and creep strain accumu- 
lations follow the same computational scheme. Sample calculated results for 
the first case analyzed are shown In table I. 

The creep analysis approach described above was Incorporated In the com- 
puter code previously developed to automate the simplified analytical proce- 
dure. This procedure Initially follows the elastic effective stress-equivalent 
strain Input until the occurrence of Initial yielding. The stress-strain solu- 
tion then proceeds along the Initial yield surface as determined from the 
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stress-strain properties. At each load step during yielding the stress shift 
(difference between new yield stress and stress predicted from elastic analy- 
sis) from the original Input data Is calculated. Elastic load reversal Is 
signaled when the Input stress Is less than the yield stress from the previous 
load step. During elastic unloading, the stresses are translated from the 
original elastic analysis solution by the amount of the calculated stress 
shift. Reverse yielding occurs when the stress reaches the reverse yield sur- 
face as determined from the hardening model. Again, the solution follows the 
yield surface until another load reversal Is Indicated when the stress based 
on the shifted elastic solution Is less than the yield stress. The elastic 
response during load reversal Is obtained by translating the original elastic 
solution according to the new stress shift calculated during reversed yielding. 
The stress-strain response for subsequent cycles Is computed by repeating this 
procedure of Identifying load reversals, tracking reverse yield surfaces and 
translating the original elastic solution during elastic loading and unloading. 
The creep computations previously discussed are performed for all load steps 
Involving dwell times. 

Neuber-type corrections for strain redistribution during mechanical load- 
ing cycles were Implemented In the plasticity calculations. The strain redis- 
tribution correction Is applied to the Ideal local total strain obtained from 
the elastic solution. These corrections were Incorporated In a version of the 
computer code using a kinematic hardening model to characterize the yield sur- 
face under cycling. A representation of a cyclic stress-strain curve by a 
bilinear kinematic hardening model Is Illustrated In figure 3. The loci of 
the tips of the cyclic curves Is described by the relation 

« = K( Cp )" (13) 


The work hardening slope, Ep, for the kinematic hardening model was deter- 
mined from energy considerations to give the same strain energy, as Indicated 
by the enclosed area In figure 3, as the actual stress-strain curve. This 
work hardening slope Is defined by 



(14) 


The code will automatically avoid these corrections for thermal loading 
problems where there are no applied mechanical loads. Without applied loads, 
the Neuber method would be Inapplicable since the stress/straln concentration 
equations would have zero net stresses and strains In the denominators. Pro- 
vision Is also made for the user to circumvent the corrections for other situ- 
ations where they would not be appropriate. These situations Include locally 
strain controlled problems and problems where the total strain Input Is based 
on strain measurements rather than elastic finite-element analyses. 


Since the stable cyclic stress-strain hysteresis loop Is a function of 
the plastic strain range. It Is necessary to Iterate between the maximum plas- 
tic strain assumed for the stress-strain relation of equation (14) and the 
resulting calculated maximum plastic strain. This Iterative process Is con- 
tinued until the specified and calculated maximum plastic strains are In 
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reasonable agreement, usually within three Iterations. Each Iteration causes 
some change In the size and shape of the cyclic stress-strain hysteresis loop. 
These changes, although usually small, make It difficult to directly compare 
solutions from the simplified procedure against nonlinear finite-element analy- 
ses since the cyclic stress-strain curves would be somewhat different. 

The computer program was validated by conducting simplified analyses for 
the benchmark notch problem and comparing the results to those from MARC non- 
linear analyses. The geometry of the benchmark notch specimen Is Illustrated 
In figure 4. This specimen was tested under Isothermal conditions as part of 
a program to provide controlled strain data for constitutive model verifica- 
tion. The benchmark notch test was conducted by fully-reversed mechanical 
load cycling at a constant temperature of 649 °C. A MARC analysis of this 
problem using kinematic hardening demonstrated excellent agreement with experi- 
mental data (ref. 9). Unfortunately, the dwell times In the cyclic tests were 
too short to achieve significant creep strains. Therefore, nonexperlmental 
cases Involving longer dwell times at the peak strain levels were devised to 
evaluate the simplified creep analysis method. Four variations of the bench- 
mark notch problem were analyzed; two with maximum creep strains of about 0.2 
percent on the tensile or loading part of the cycle, a third with a maximum 
creep strain under 0.1 percent on the compressive or unloading part of the 
cycle and a fourth case Involving combined tensile and compressive creep to 
Illustrate a cyclic creep problem. Cyclic yielding was described by a kine- 
matic hardening model with cyclic stress-strain data for Inconel 718 alloy 
(ref. 6). Nonlinear and elastic MARC analyses of this problem were performed 
using approximately 600 triangular plane strain elements to model a quarter 
segment bounded by planes of symmetry as shown In figure 5. The MARC solutions 
shown for the benchmark notch specimen were computed at the closest Gaussian 
Integration point to the root of the notch. 


DISCUSSION OF ANALYTICAL RESULTS 

Analytical results for the benchmark notch cases using the Improved creep 
analysis Incorporated In the simplified Inelastic procedure are discussed 
herein. Comparisons are made with MARC finite-element solutions for the same 
problems. Stress-strain cycles used for comparison purposes are In terms of* 
von Mlses effective stresses and equivalent strains with assigned tensile or 
compressive signs. The entire discussion Is based on the critical location, 
which was at the notch root. 

The conditions of the four analytical cases studied are summarized In 
table II. The analyses for the first three cases were conducted by enforcing 
the same Initial stress at the start of creep on the simplified procedure as 
was calculated for the comparable MARC analysis. This was necessary because 
small differences In the cyclic stress- strain curves could cause differences 
In the Initial stresses of about 3 percent and make It difficult to directly 
compare both creep analysis methods. 

In cases 1 and 2, 4-hr dwell times were Imposed at the maximum tensile 
load, or the maximum total strain part of the cycle. These two cases differ 
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In the load level at the start of creep. In both cases the MARC finite-element 
analyses broke down after reaching a creep microstrain of 1000 (or less than 
1 hr of accumulated creep time) because of Inability to maintain the relatively 
tight tolerance on the stress change per increment. Opening up the tolerance 
by an order of magnitude enabled the finite-element creep computations to be 
extended to 1.5 hr of accumulated dwell time. However, the looser tolerance 
resulted In a more erratic calculated stress relaxation response until the 
finite-element creep analysis was automatically terminated due to Instability. 
In contrast, the simplified analysis exhibited no computational problems and 
gave a stable solution for the full 4-hr dwell times. Comparison of the sim- 
plified creep analyses with the MARC analyses, as far as they could be carried 
out, are shown In figures 6 (a) and (b) for cases 1 and 2 respectively. Excel- 
lent agreement Is seen between the two analytical methods for both cases. 

In case 3, a 4-hr dwell time was applied at the minimum load In the cycle. 
Since the stress level at the start of creep was much smaller than In the first 
two cases, there was no computational problem with the finite-element creep 
analysis for case 3. Again, excellent agreement was shown between the simpli- 
fied and finite-element creep analyses as Is demonstrated In figure 6(c). In 
all three cases, the simplified solution agreed with the finite element results 
within 1 percent In terms of stress during relaxation and 6 percent In terms 
of creep strain. 

The ability of the simplified method to calculate cyclic creep was evalu- 
ated In case 4. Dwell times were applied at both the maximum and minimum 
cyclic load levels. The dwell time at the maximum load was limited to 1 hr to 
avoid the computational Instability problem with the finite element analysis 
at longer dwell times. A 10-hr dwell time, rather than 4 hr as In case 3, was 
applied at the minimum load to obtain larger creep strains In the minimum 
strain part of the cycle. A comparison of the simplified and finite-element 
solutions for case 4 Is shown In figure 7. The comparison was limited to the 
first cycle because of the difficulty of maintaining computational stability 
during the second cycle for the finite element analysis. Calculated points 
from the simplified analysis are shown only for the beginning and end of each 
dwell time. There Is a noticeable difference between the creep behavior during 
the compressive dwell time shown In figure 7 and that Illustrated In figure 1. 
This difference Is due to the fact that Since the finite-element analysis con- 
siders the equivalent creep strain change during compression as a positive 
number, the magnitude of the total strain during the second dwell period can 
only be reduced. To maintain consistency with the MARC creep analysis, the 
simplified analysis was adapted to have the same type of creep behavior on 
load reversal . However, for an actual structural analysis of a cyclic creep 
problem the simplified procedure would consider a creep strain change during 
compression as a negative number and give the behavior shown during the com- 
pressive dwell time In figure 1. 

Most of the discrepancies between the simplified and finite-element hys- 
teresis loops In figure 7 were due primarily to the Iteration process built 
Into the simplified analysis computer code which resulted In the use of a 
cyclic stress-strain curve slightly different from that used in the MARC 
elastic-plastic analysis. During the Initial monotonic loading, the two 
analytical methods agreed, as In cases 1 to 3, within 1 percent in relaxed 
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stress and 6 percent In creep strain. However, during load reversal dwell 
time the differences Increased to 1.8 percent In terms of relaxed stress and 
almost 2 to 1 In terms of cumulative creep strain. The larger disagreement In 
creep strain was mainly due to the accumulation of error as one proceeds 
through the cycle and the sensitivity of creep calculations to relatively small 
discrepancies In the stress levels. The simplified analysis of the benchmark 
notch problem for case 4 used 0.3 percent of the central processor unit (CPU) 
time required for the MARC nonlinear analysis of the finite element model shown 
In figure 5. The MARC computing time could have been reduced by concentrating 
on the notch region and applying proper boundary conditions. However, the 
simplified procedure would still be much faster because It focusses on a single 
point. 


SUMMARY OF RESULTS 

An Improved uniaxial creep analysis method was developed for use with a 
simplified Inelastic analysis procedure for calculating the stress-strain 
history at the critical location of a thermomechanical ly cycled structure. 

This creep method used a self-adaptive time Incrementing scheme for subdividing 
the dwell time period Into smaller Increments for creep analysis. In this 
approach, cumulative creep was the sum of the Initial creep strain, the recov- 
ery from stress relaxation and the Incremental creep strain. Neuber-type cor- 
rections were applied to account for strain redistribution and residual 
stresses due to plastic strain reversals. Analytical predictions from the 
simplified procedure for benchmark notched plate problems were compared to 
nonlinear finite-element solutions. The following general conclusions were 
drawn from the evaluation 
of the method: 

1. The predicted stress-strain response using the simplified procedure 

with the uniaxial creep analysis showed excellent agreement with three- 
dimensional elastic-plastic-creep finite element solutions using the MARC 
program. In all the monotonic creep analysis cases studied, the simplified 
and finite-element solution for the critical crack Initiation location agreed 
within one percent In terms of relaxed stress and 6 percent In terms of 
cumulative creep strain. * 

2. The simplified creep analysis was extremely stable. Creep calculations 
could be carried out to tight tolerances on stress changes and continued under 
conditions where the nonlinear finite-element analysis became computationally 
unstable and broke down. 

3. The simplified procedure computed stress-strain hysteresis loops for 
the critical location using 0.3 percent of the CPU time required for elastlc- 
plastlc-creep finite-element analyses of the full finite element model. 
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APPENDIX A 

DERIVATION OF BASIC CREEP EQUATIONS 

Equation (2) Is obtained by combining equations (1), (10), (11) and (12). 
From equatlon(lO) 

e c( 1+1 ) = C c 1 + Ac c 1 - e 1 c cR1 

Substituting equation (1) for e c ^, equation (11) for Ae c ^ and 
equation (12) for g c r^ yields 


C c( 1 +1 ) = 



+ 




C 

1+1 
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e^(c 


B 

1=1 



c 
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or 

'c(i.i) - (a) B {<’ ‘li } 

From figure 2, the creep strain at the end of an Increment Is given by 
B C B C 

C c(1+1) = (a) tf i+i [t i+i + Y i A V + (a) °1+1 [t 1 + (1 + V A V 
Equating the above two expressions leads to the form of equation (2) 



Other equations given In the Analytical Procedure can be obtained from 
similar mathematical manipulations. 
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TABLE I. - SAMPLE CALCULATION (CASE 1) 


Inc 

ti > 
hr 

At-j , 
hr 

t-j+1, 

hr 

Ae ci 

8 i e cRi 

£ ci >i 
xlO 4 

Total equivalent strain, 
E c( i+1) 

Effective stress, 
°i+l,MPa 

Simplified 

MARC 

Simpl i f i ed 

MARC 

1 

0 

0.002767 

0.002767 

1.385xl0' 4 

0 

0 

1.385xl0 -4 

1.402x10-4 

809.7 

809.0 

2 

0.002767 

.002767 

.005535 

.478 

0.0444 

1.385 

1.818 

1.793 

807.8 

810.1 

3 

.005535 

.004151 

.009686 

.501 

.0213 

1.818 

2.298 

2.323 

806.1 

807.1 

4 

.009686 

.006226 

.01591 

.554 

.0253 

2.298 

2.827 

2.826 

804.1 

804.5 

5 

.01591 

.00934 

.02525 

.6262 

.0358 

2.827 

3.417 

3.396 

802.0 

801.7 

6 

.02525 

.01401 

.03926 

.7201 

.0430 

3.417 

4.086 

4.031 

799.4 

798.7 

7 

.03926 

.02101 

.06028 

.8292 

.0745 

4.086 

4.841 

4.751 

796.4 

795.4 

8 

.06028 

.03152 

.09180 

.9581 

.1464 

4.841 

5.652 

5.565 

792.3 

791.9 

y 

.09180 

.04728 

.1391 

1.098 

.1917 

5.652 

6.566 

6.490 

787.4 

788.1 

10 

.1391 

.0709 

.2100 

1.246 

.256 

6.566 

7.555 

7.536 

783.6 

784.0 

11 

.2100 

.1064 

.3164 

1.447 

.243 

7.555 

8.759 

8.723 

779.7 

779.7 

12 

.3164 

.1596 

.4760 

1.669 

.301 

8.759 

10.127 

10.060 

775.5 

775.2 

13 

.4760 

.2394 

.7153 

1.926 

.361 

10.127 

11.692 

11.580 

771.4 

770.6 

14 

.7153 

.3591 

1.0744 

2.199 

.541 

11.692 

13.352 


766.1 

770.6 

15 

1.0744 

.5387 

1.6131 

2.527 

.532 

13.352 

15.347 


761.7 


16 

1.6131 

.8081 

2.4212 

2.905 

.625 

15.347 

17.627 


757.1 


17 

2.4212 

1.2122 

3.6334 

3.334 

.734 

17.627 

20.228 


752.5 


18 

3.6334 

.3666 

4.0000 

.845 

.746 

20.228 

20.327 


752.1 



at^ = 0.05 hr 
t = 4 hr 
a = 0.01 


TABLE II. - ANALYTICAL CASES 


Case 

Load applied during 
dwell time, 
kN 

Owel 1 
time, 
hr 

Tolerance coefficient 

on stress change 

Simplified analysis 

MARC analysis 

1 

33.3 

4 

0.01 

0.01 

2 

35.1 

4 

.01 

.01 

3 

-17.3 

4 

.01 

.01 

4 

' 35.1 

1 

.01 

.10 


-17.3 

10 

.01 

.10 




CREEP 



STARTS 

Figure 1. - Typical stress-strain cycle with creep dwell times. 



Figure 2. - Modified time hardening law. 




Figure 4. - Benchmark notch specimen (x t = 1.9). (All dimen- . 
sionsinmm. ) 
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